Vortex precession in Bose-Einstein condensates: observations with filled and empty 

cores 
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We have observed and characterized the dynamics of singly quantized vortices in dilute-gas Bose- 
Einstein condensates. Our condensates are produced in a superposition of two internal states of 
87 Rb, with one state supporting a vortex and the other filling the vortex core. Subsequently, the 
state filling the core can be partially or completely removed, reducing the radius of the core by as 
much as a factor of 13, all the way down to its bare value. The corresponding superfluid rotation 
rates, evaluated at the core radius, vary by a factor of 150, but the precession frequency of the 
vortex core about the condensate axis changes by only a factor of two. 

PACS number(s): 03.75.Fi, 67.90,+z, 67.57.Fg, 32.80.Pj 



The dynamics of quantized vortices in superfluid he- 
lium and superconductors have been fascinating and im- 
portant research areas in low-temperature physics |Q. 
Continued study of vortex dynamics may, for example, 
lead to a better understanding of energy dissipation in 
these systems g|. Work on optical vortices has also 
become an active area of research ||. More recently, 
demonstrations of the creation of quantized vortices in 
dilute-gas Bose-Einstein condensates (BEC) ||,[| have 
proven to be striking examples of the similarities between 
the condensed matter, optical, and dilute-gas quantum 
systems. Because of the observational capabilities of 
dilute-gas BEC experiments and the ability to manipu- 
late the quantum wavefunction of the condensates, these 
systems provide a unique approach to the study of quan- 
tized vortices and their dynamics. This paper reports di- 
rect observations and measurements of singly quantized 
vortex core precession in a BEC. 

Numerous theoretical papers have explored the ex- 
pected stability and behavior of vortices in BEC [p|~[i3| . 
One interesting predicted effect is vortex core precession 
about the condensate axis P,p|-|l2|| . Radial motion of the 
core within the condensate can also occur, and may be 
understood as being due to energy dissipation and damp- 
ing processes. 

Core precession may be described in terms of a Magnus 
effect - a familiar concept in fluid dynamics and super- 
fluidity jlj]. An applied force on a rotating cylinder in a 
fluid leads to cylinder drift (due to pressure imbalances 
at the cylinder surface) that is orthogonal to the force. 
Analogously, a net force on a vortex core in a superfluid 
results in core motion perpendicular to both the vortex 
quantization axis and the force. In the condensate vortex 
case, these forces can be due to density gradients within 
the condensate, for example, or the drag due to ther- 
mal atoms. The density-gradient force may be thought 
of as one component of an effective buoyancy: just as 
a bubble in a fluid feels a force anti-parallel to the lo- 



cal pressure gradient, a vortex core in a condensate will 
feel a force towards lower condensate densities. The to- 
tal effective buoyancy, however, is due less to displaced 
mass (the "bubble" ) than it is to dynamical effects of the 
velocity-field asymmetry, which in turn is a consequence 
of a radially offset core. Typically, the total buoyancy 
force is towards the condensate surface, and the net ef- 
fect is a precession of the core around the condensate axis 
via the Magnus effect. Drag due to the motionless (on 
average) thermal atoms opposes core precession, causing 
the core to spiral outwards towards the condensate sur- 
face. In the absence of this drag (for temperature ~ 0), 
radial drift of the core may be negligible. 

Our techniques for creating and imaging a vor- 
tex in a coupled two-component condensate are de- 
scribed in refs. ||,|l4|. The two components are the 
\F = l,mp = —1) and \F = 2,mp = 1) internal states 
of the 87 Rb atom, henceforth labeled as states |1) and 
|2) respectively. We start with a condensate of 10 6 
1 2) atoms, confined in a spherical potential with oscil- 
lator frequency 7.8 Hz. A near-resonant microwave field 
causes some of the |2) atoms to convert to |1) atoms. The 
presence of a rotating, off-resonant laser beam spatially 
modulates the amplitude and phase of the conversion. 
The net result is a conversion of about half of the sam- 
ple into an annular ring of |1) atoms with a continuous 
quantum phase winding from to 2tt about the circum- 
ference - a singly quantized vortex. The balance of the 
sample remains in the non-rotating |2) state and fills the 
vortex core. With resonant light pressure we can selec- 
tively remove as much of the core material as we desire. 
In the limit of complete removal, we are left with a single- 
component, bare vortex state. 

In this bare-core limit, the core radius is on the order 
of the condensate healing length £ = (87rn a) -1 / 2 , where 
uq is the peak condensate density and a is the scattering 
length. For our conditions, £ = 0.65 /im, well under our 
imaging resolution limit. The bare core can be observed 
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after ballistic expansion @] of the condensate, but this is a 
destructive measurement. On the other hand, if we leave 
some of the |2)-state atoms filling the core, the pressure 
of the filling material opens up the radius of the |1) vortex 
core to the point where we can resolve the core in a time 
series of nondestructive phase-contrast images. 

Filled- core dynamics - We first discuss vortex dynam- 
ics in two-component condensates, where 10% — 50% 
of the atoms were in the |2) fluid filling the |1) vortex 
core. We took successive images of the |1) atoms in 
the magnetic trap, with up to 10 images of each vor- 
tex. The vortex core is visible as a dark spot in a bright 
|1) distribution, as shown in Fig. 0(a). Instabilities in our 
vortex creation process usually resulted in the creation of 
off-center vortex cores, allowing us to observe precession 
of the cores. We observed precession out to ~ 2 s, af- 
ter which the |2) fluid had decayed to the point that the 
vortex core was too small to be observed in the trapped 
condensate. 

The recorded profile of each trapped condensate was 
fit with a smooth Thomas-Fermi distribution. Each vor- 
tex core profile was fit with a Gaussian distribution to 
determine its radius and position within the condensate. 
From the fits, we determined the overall radius Rt of the 
trapped condensate (typically 22 /um), the HWHM ra- 
dius r of the filled vortex core, and the displacement dt 
and angle 9 t of the core center with respect to the con- 
densate center. Core angles and radii for the images in 
Fig. ^(a) are shown in Figs. |l|(c) and (d). The vortex core 
is seen precessing in a clockwise direction, which is the 
same direction as the vortex fluid flow around the core. 
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FIG. 1. (a) Seven successive images of a condensate with 
a vortex and (b) their corresponding fits. The 75-/xm-square 
nondestructive images were taken at the times listed, refer- 
enced to the first image. The vortex core is visible as the dark 
region within the bright condensate image, (c) The azimuthal 
angle of the core is determined for each image, and is plotted 
vs. time held in the trap. A linear fit to the data indicates a 
precession frequency of 1.3(1) Hz for this data set. (d) Core 
radius r in units of healing length £. The line shown is a linear 
fit to the data. 



The angular precession frequency was determined from 
the time dependence of 9 t (Fig. |l|(c)). This and other sim- 
ilar data sets showed no reproducible radial motion of the 
core over the times and parameters examined. However, 
consistent decrease in the size of the core was observed, 
which we interpret as being due to known decay of the 
2) fluid through inelastic atomic collisional processes. 

For each data set, we determined a mean core radius 
and displacement. The data cover a range of core radii 
(r = 7£ to 130, displacements (dt = 0.l7R t to 0.48i? t ), 
and percentage of atoms in the core (10% to 50%). Ex- 
cept for a few "rogue vortices" (discussed below), the 
measured precession frequencies are clustered around 1.4 
Hz, as shown in Fig. ^, precessing in the same direction 
as the vortex fluid rotation. The data (Fig. 0(a)) sug- 
gest a slight increase in frequency for cores further from 
the condensate center, in qualitative agreement with two- 
dimensional numerical simulations for two-component 
condensates |l5[. We also see (Fig. ||(b)) a slight decrease 
in precession frequency for larger cores. 

As indicated in Fig. ||, a few vortex cores exhibited 
precession opposite to that of the fluid flow, but with ab- 
solute frequencies similar to the bulk of the data. The 
quality of the corresponding vortex images was routinely 
lower than for the positive-frequency precession points, 
with vortices looking more like crescents and "D" shaped 
objects rather than like the images of Fig. § We spec- 
ulate that this "inverse precession" may be due to dis- 
torted configurations of vortices in the two-component 
condensate in which vortex dynamics are not as well un- 
derstood. Such configurations have received recent theo- 
retical attention Jl3] ], 
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FIG. 2. Compiled data for filled vortex core precession, 
with each data point extracted from a series (as in Fig. |l]) 
of nondestructive images of a single vortex. Precession fre- 
quency is plotted vs. (a) core displacement d t in units of con- 
densate radius R t , and (b) core radius r in units of healing 
length £. Circles correspond to positive frequencies and filled 
squares to negative frequencies. (Positive frequency is defined 
as core precession having the same handedness as the vortex 
angular momentum.) The triangle at r = £ shows for refer- 
ence the average measured precession frequency of many bare 
vortices (see text and Fig. 3|(b)). A line is drawn as a guide 
to the trend in frequency vs. core size. 
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Bare core dynamics - To examine the dynamics of 
bare vortices, our procedure consisted of taking a non- 
destructive phase-contrast picture of the partially filled 
|1) vortex distribution (Fig. ||(a) inset), as previously dis- 
cussed, followed by complete removal of the core filling 
Jl6| . We then held the bare vortex in the trap for a vari- 
able hold time th, after which the condensate was released 
from the trap. We took a final near-resonance phase- 
contrast image 17| of the atomic distribution (Fig. ||(a)) 
after the condensate had ballistically expanded by a fac- 
tor of ~ 3.5 jl8| and the core had expanded jl9| to a fit 
radius of ~ 9/im. 

Displacements dt and angular positions 9t of the cores 
for the in-trap images were extracted as described before. 
The images of the expanded clouds were fit with identical 
distributions, and the Thomas-Fermi radius R e of the 
expanded cloud and the vortex core displacement d e and 
angle 9 e were obtained for each image. For each pair 
of images, we determined the angular difference A9 et = 
8 e — 6 t between the cores in the expanded and in-trap 
images. We also determined the core displacement ratio 
d e /dt, an indicator of the radial motion of the core during 
the hold time th- 

From the measurements of A9 et at different hold times 
th (Fig. |^(b)), we find a bare core precession frequency 
of 1.8(1) Hz, which is slightly faster than the precession 
of filled cores and consistent with the trend shown in 
Fig. 0(b) for filled cores. To emphasize that our mea- 
surements of filled and empty cores are different limits in 
a continuum of filling material, we indicate the measured 
bare core precession frequency in Fig. ^|(b) as a point at 
r = £. 

From Fig. |^(b) it is apparent that the structure and 
content of the vortex core have a relatively modest effect 
on precession frequency. One can calculate, for instance, 
the fluid rotation rate v r at the inner core radius. The 
value of v r is given by the quantized azimuthal superfluid 
velocity evaluated at the radius of the core, divided by the 
circumferential length at that radius. For the bare-core 
vortex (the triangle in Fig. ||(b)) v r is about 260 Hz, while 
for the largest filled cores of Fig. ^(b) (for which nearly 
half of the sample mass is composed of core filling) , v r is 
only about 1.7 Hz. Thus between vortices whose inner- 
radius fluid rotation rates vary by a factor of 150, we see 
only a factor of 2 difference in precession frequency. 

The slower precession of filled cores can be understood 
in terms of our buoyancy picture. Due to its slightly 
smaller scattering length, |2) fluid has negative buoyancy 
with respect to |1) fluid, and consequently tends to sink 
inward towards the center of the condensate [^0|. With 
increasing amounts of |2) material in the core, the inward 
force on the core begins to counteract the outward buoy- 
ancy of the vortex velocity field, resulting in a reduced 
precession velocity. It is predicted that with a filling ma- 
terial of sufficiently negative buoyancy in the core, the 
core precession may stop or even precess in a direction 
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FIG. 3. (a) Ballistic expansion image of a vortex after all 
1 2) atoms have been removed. The dark spot is the bare vor- 
tex core. Inset - the corresponding, preceding in-trap non- 
destructive image of the partially-filled core, (b) Angular 
differences A# e i between vortex cores from the in-trap and 
expansion images, plotted against hold time th in the mag- 
netic trap. The line is a fit through the data, indicating a bare 
core precession frequency of 1.8(1) Hz. (c) Radial core motion 
is determined by (d e / R e ) / (dt / Rt) , the ratio of the fractional 
core displacements from the expansion and in-trap images of 
each data set. The data are shown as open circles, with the 
average of all data at each given hold time plotted as a filled 
triangle, (d) Survival probability of a vortex, defined as the 
conditional probability for observing a vortex in an expanded 
image given the observation of a vortex in the corresponding, 
pre-expansion in-trap image. Bare cores are not observed for 
hold times th > 1.1 s. 

opposite to the direction of the fluid flow jl5|], but our 
data do not reach this regime. 

Various theoretical techniques involving two- and 
three-dimensional numerical and analytical analyses have 
been explored to calculate the precession frequency of a 
vortex core within a condensate ]Tc| ^T^j2T[ | . We briefly 
compare those most readily applied to our physical pa- 
rameters, assuming a spherical, single-component con- 
densate with 3 x 10 5 atoms in a non-rotating trap. Where 
relevant, we assume a core displacement of d t = 0.35i?t 
which is near the center of our measurement range. 
A two-dimensional hydrodynamic image charge analysis 
has been analytically explored in the non-interacting |l2j ] 
and weakly-interacting (Thomas-Fermi) pl[ limits. The 
latter of these predicts a bare core precession frequency 
of ~ 0.8 Hz. Svidzinsky and Fetter's two-dimensional 
|llf and three-dimensional [0] solution to the Gross- 
Pitaevskii equation predicts a precession frequency of 
~ 1.3 Hz. Jackson et. at jl0| have obtained results in 
close agreement with this analytical solution using a nu- 
merical solution to the Gross-Pitaevskii equation. Fi- 
nally, a two-dimensional simulation by McGee and Hol- 
land p5| using a steepest-descents technique predicts a 
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precession frequency of 1.2 Hz. 

Measurement of d e /dt for different hold times th shows 
the radial motion of the bare cores and is a probe of en- 
ergy dissipation of the vortex states. The plot of Fig. ||(c) 
displays no trend of the core towards the condensate sur- 
face during th, indicating that thermal damping is neg- 
ligible on the 1 s timescale p^ |. However, a peculiar 
feature of our bare core measurements is the lack of visi- 
ble cores in the expanded pictures for hold times of th > 
1 s (Fig. ||(d)). The absence of observed radial motion of 
the core suggests that the vanishing of the vortex is not 
due to thermal damping (which would result in the core 
moving towards and annihilating at the condensate sur- 
face) . Possible explanations are that the vortex core may 
tilt or bend along its length Q , suppressing contrast in 
optical depth below our signal-to-noise threshold. 

The hypothesis that the vortices are disappearing due 
to imaging limitations rather than true decay is sup- 
ported by an alternative assay for the presence of vor- 
tices: It has been predicted that the presence of a vortex 
will break the degeneracy of the m = ±2 condensate sur- 
face modes [p4[ . We have experimentally verified this 
prediction p5| and used this technique to verify that the 
vortex continues to be present in the condensate well af- 
ter we can no longer image a core in expansion. These 
measurements will be the topic of a future paper pcfl . 

Through a combination of destructive and nondestruc- 
tive imaging techniques we have obtained measurements 
of vortex dynamics in bare- and filled-core vortices in 
dilute-gas BEC. Vortex precession frequencies show only 
modest dependence on the radius and content of the vor- 
tex core. We plan to extend these studies to investigate 
higher-order dynamical behavior and to characterize the 
dissipative effects of finite temperatures. 
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ray Holland and Sarah McGee. This work was supported 
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